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Abstract 

An improved approximation via Havriliak-Negami (HN) functions to 
the Fourier Transform (FT) of certain Weibull distributions, —'i/’/s, (the 
time derivative of the Kohlrausch-Williams-Watts function), is given for 
a large interval of frequencies: aj/27r £ [0,10^^] if 0 < /3 < 1 a-nd aj/27r £ 
[0,10^] if 1 < P < 2. The model is free from the usual numerical 
distortions, or restrictions associated to sampling step and hnite size, 
present in similar adjustments to complex relaxation functions. Fur¬ 
ther indicates that the identification of (FT) Weibull data with a dou¬ 
ble HN approximant, -i/i/j ~ Ap 2 HN, is quite exact locally even though 
the parameters involved should vary adiabatically with the frequency, i.e. 
{ai, 2 , 7 i, 2 ,'ri, 2 , A}(cu). This fact is the base for the high sensibility of the 
parameter models, as functions of /3, to the available data and sampling 
associated errors. Presumably is also the root of the different propos¬ 
als for asymptotic traits of a ■ '){P)i, 2 \lj^oo functions found in scientihc 
literature. 


Introduction 

The Kohlrausch relaxation function, 4>K,i3it) = exp—0 < t < oo, /3 < iQ 
initially proposed to explain the discharge of capacitors [1], has found later 
several applications in many fields not directly related to electrostatics. And 
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while it was half-forgotten in the scientific literature which was more worried 
about practical applications, it did not dissapear at any time of the mathe¬ 
matical heritage due to the deep connections that such function and its se¬ 
quel, /3 > 1, have with number theory m El a 0 El [T], stable distributions in 
probability n a uni m EH n HI, theta function extensions ESI EH, or the 
special functions obtained from certain kernels of Mellin-Barnes type transfor¬ 
mations employed in fractional calculus EH EH EH IMl mi HI However 
it is a result of its use by Williams and Watts (23], in the field of character¬ 
ization of dielectric properties of polymers |24l l25l [26], that its employment 
became ubiquitous in many areas of Physics and Chemistry. Hereafter not only 
associated with specific electrical phenomena but a plethora of them in the 
most diverse fields. It emerges in luminiscence [23 [251 [221 EH Ell EH ESj) su¬ 
percooled glasses EH EH ESI) capacitors El EH EH EH, rheology EHSHIIH, 
spin glasses [34], torsion of galvanometric threads [IH, dielectric spectroscopy 
(^ EH im SH BH113 SHI SH, magnetic resonance (501 EH EH, autocorrelation 
functions in molecular dynamics [23 EH EH EH, econophysics |54[ EH , protein 
folding [2HEHEHE3EH] and even astrophysics |22| - 

Since it shows up in many branches of physics of complex systems |38l 1581 
EH EH EH and soft matter [a SH S3 SH 0 EH EH , among others, seems not 
be linked to a precise interaction or a particular physico-chemical phenomenon, 
instead it is more appropriate to associate it to an emergent property. In such a 
case any ensemble of interacting elements organized in multiscale clusters whose 
local relaxations, or restructuring bonds, proceed to jumps with random wait 
times of type t~°‘ should present an autocorrelation, or decay, of Kohlrausch’s 

nature |3H SH SH EH EH E3 EH EH- Thus its hierarchical structure, and 

possibly experimental difficulties, advises us to study this relaxation not only 
through its response along time but also in other spaces of representation, such 
as of frequency [33 SH EH EH EH EH EH EH EH ■ An exhaustive account of the 
distribution of intermolecular configurational energies, characteristic relaxation 
times -real or virtual-, and frequency response will give us a valuable set of 
techniques to tie in different analytical and laboratory procedures. And indeed 
such possibilities are available since the Transforms of Fourier and Laplace, 
(w—space and r—space respectively), exist for the case 0 < /3 < 2 [51 121 [TOl [T21 
EH- Either as a convergent series or an asymptotic one, according to subcases 
0 < /3 < 1 and 1 < < 2 and depending on we are talking about high or low 

frequency [HIHEHEHEHEH- 

Unfortunately these series present several problems of convergence mainly 
near zero, and in particular the best known of them shows up an essential 
singularity and the others are simply asymptotic dH EH EH EH EH EH , i. e. non 
convergent as a whole, and with coefficients associated to the special function 
gamma, F, over the entire reals which makes calculations more tedious than in 
other simpler expressions [31EHEHE3- This makes the sum of such series 
a burdensome task which is possible to get round with numerical methods or 
resummation of series [13 EH EH EH- Nowadays both techniques allow a direct 
and quick calculation of Laplace and Fourier Transforms with enough accuracy 
to be of great utility in spectral filtering, or analysis, of laboratory data [IHEH]. 
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Nevertheless could be missed a concise mathematical formula to give account, 
even approximately, of the mentioned transforms to accelerate the calculations 
or evaluate repeatedly such functions [53j [65j |7D]. The latter is usually a must in 
optimization since the multiple evaluation of objective functions with variational 
parameters increases drastically computational loads. Genetic algorithms are a 
good example of this problem. On the other hand a short analytical expression 
makes it easier to compare with the most common mathematical functions in 
the complex domain which are already used in the field of dielectric relaxation, 
such as Debye [80], Cole-Cole |8T], Cole-Davidson [82], Havriliak-Negami |88l[84] 
and others. 

We use the notation xpi^) = /g for the one-sided Fourier 

transform of the Kohlrausch function and '4’pi!^) = ~Jj(/>ic,/ 3 (t)dt for 
minus the transformation of the Weibull distribution. They are both related 
in w—space by the algebraic closure relationship "(ppioj) iuJX 0 (i^) = 1. In this 
way it is indistinct to work with ijjp or Xi 3 i Ih® former is however preferable 
since it allows to remove a pole and obliterates any singularity in the neigh¬ 
borhood of a; « 0 while using the resummation formula for the original series 
of inverse powers of which describes x^, (for more details see Ref. |SI|). 
Other advantages of employing ipp are homogeneity, boundedness and similar¬ 
ity of this family of functions if the shape parameter j3 is varied, as we pointed 
out previously m- We will improve this reference, and others similar in nature 
isaiMjiiniini, extending the limits imposed by numerical samplings in t—space 
and making consistent the description of high frequencies decays for ipp with 
a unique set of strict Havriliak-Negami functions when /3 < 1, and with set of 
same kind although parametrically extended if /3 > 1 [^ . Besides we will show 
how the apparently supernumerary parameters of the Havriliak-Negami func¬ 
tion |88l IMj . ; 0 < a ,7 < 1, when approximating the mentioned 

transform are uniquely determined by the parameter /3. Thus these models with 
such strong dependency {a, 7 ,r}(/?) are equated with others quite new of less 
parameters |65| . 


1 Analytical considerations 

1.1 The asymptotes of the data 

As it was shown |53| . a double approximant of Havriliak-Negami functions de¬ 
scribe fairly well the FFT of the derivative of the Kohlrausch function as well as 
the Cole-Davidson-Kohlrausch family |5^ I65| . Being the main source of error 
an small impairment of the fit in the low frequencies. Nevertheless as the FFT 
is obtained from the original function by means of a finite window sampling, 
which in the Fourier space is reflected as a convolution with a sine function, 
a round off of the spectra in the high frequency border of the domain is ob¬ 
tained. Consequently the asymptotic behaviour of the original function and 
the approximant will differ, at least in the high frequency domain limits. Be¬ 
sides, due to the smearing of the peaks when the numerical transformation is 
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involved, a difference between the theoretical FT and that of the numerically 
sampled function is expected. This is why the double approximant fitted to the 
numerical transformation possibly increments even more than its mathemati¬ 
cal limitations the difference with the FT of at the low frequencies. 

The question raised is then how sensitive are the parameters obtained in the 
optimization to this deformation, inasmuch as the asymptotic behaviour of the 
Kohlrausch’s (or Weibull’s [551186] 1 Fourier Transform has been indexed to them 
{v.gr. ai ■ ji = ai- 7i(/3) and the rest are functions of /3 too) [751 ITTl 157] . 

To assess the influence of finite-size deformation in the parameters obtained 
by fitting a couple of Havriliak-Negami functions to the Fourier Transform, 
ippluj), we have performed two different simulations of it with a general purpose 
mathematics package (Mathematica"^^). The aim is to obtain accuracy and 
avoid numerical oscillations as far as possible, the latter a characteristic of the 
integrals involved that makes difficult results at very high frequencies, as well 
as to get the proper asymptotic behaviour of the tails. 

As it was for the FFT case the first series comprehends frequencies from 
u = 0 to v = 500.0005 in steps of 6 iy = 1/999.999, being w = 27ri/ [53]. The 
range of parameter /3 simulated is (0, 2.00], and two subintervals of different 
traits are analyzed independently. One /3 s(0,1.00] the other /? s(1.00,2.00] with 
test points chosen as {Q,l}.xy with x = {0,1,2,3,4,5,6,7,8,9} and y = (0,2,5,8} 
and the end point (3 = 2.00. 

For the second series the /3 points are the same aforementioned but the choice 
for the interval of frequencies is different for each of the cases /3 <1 and /3 >1. 
So if P <1, V g[ 0,10^^] and for /3 >1, v g[ 0,10^] and the reason for these distinct 
intervals is the increasing numerical noise that overshadows the signal. This 
makes any calculation useless beyond u = 10^ for /3 >1 and mostly for /3 >1.90. 
Besides the steps of frequencies are not homogeneous in the whole interval as 
they were in the first series; they are now in this way only in logarithmic scale. 
This is, in the interval v g[ 0,10“®] the step takes the value 5v = 10"®, and in each 
interval v €(10“b 10“’“'"^] the increment is 5ui = 9 * 10“*“® with ai any integer 
of set (-5,-4,...,-1,0,1,...,11,12} if /? <1, or set (-5,-4,...,-1,0,1,...,6,7} if /3 >1. So 
with this second series it is possible to sample a wider domain not affordable 
with the fine grain step we used in the first series, the drawback here is the 
different implicit ’weight’ the points acquired in both series when the fitting to 
a double Havriliak-Negami approximant is made. 

The reason for such a large domain of frequencies it is not the reconstruc¬ 
tion purpose, as was shown that a domain of low to medium frequencies, (i. e. 
V g[ 0,500.0005]), were enough to depict the function e“* quite accurately by 
means of inverse FFT of xpi'^) E3]- Nor, of course, the physical reconnaissance 
by means of dielectric spectroscopy as usually is done in polymer or metallurgi¬ 
cal sciences because just one type of relaxation is rarely presented in isolation. 
Usually different kind of relaxations, associated to several size scales and diverse 
phenomena, are registered together in contiguous domains of frequency spoiling 
each diagram the head and tail of its neighbors. Rather we are interested in 
tails from a mathematical point of view, since a proper knowledge of the func- 


4 


tion -and this includes the exact tail shape-, will allow us to concoct the more 
precise approximant possible in the whole domain, and which should include 
the difficult-by-series description at low frequencies, i.t. near u ^ Q. 

The modulus of function (w) presents the following asymptotic behaviour 
in the domain of frequencies: when w —>■ 0 and |' 0 ,g(a;)|~r(/? -|- l)/a;^ 

as w —t oo being monotonically decreasing in the intermediate values. Never¬ 
theless the interval of frequencies in which a significant amount of this drop¬ 
ping takes place depends strongly with the value of parameter /3. While for 
values near to zero, /3 —> 0, the decrement of the signal is slow and mild, 
near to one, /3 ~ 1, is quicker and sharper. For example in the interval 
V S [0.001,500.0005], the initial values of |V'/ 3 (w)| are: (0.665237, 0.793929, 
0.999607, 0.999973, 0.999980) and the respective percentage drops are (14.462%, 
56.382%, 98.429%, 99.931%, 99.968%), at the end of such domain, when /? is 
set to (0.02, 0.10, 0.50, 0.90, 1.00). This suggests important features about the 
functions \il}p{(jj)\, namely the plateau near z/ ~ 0 is growing in length as /? —t 1 , 
the mentioned evolution towards a steeper function, and consequently the less 
important contribution of the tails to the decrement as P varies. Such behaviour 
is held and even exacerbated in the interval 1 < /3 < 2 , not only because the 
drop is greater-(99.9851%, 99.9992%, 99.9999%) for/? e{1.10, 1.50, 2.00}, with 
starting values (0.999985, 0.999992, 0.999996) of |'0^(a;)|- but also a sudden 
change in the decreasing of the functions is perceived. 



Figure 1: a) Logarithmic derivative of logjg 11 / 1 , 3 ( 0 ;)| vs. logj^g(o;/27r). b) Logarithmic second 
derivative of logj^p |i/;^(a;)| vs. logjg(o;/27r). c) Exponent a of fitting |^^(a;)| ~ 0 ;““ to data 
tail. Subindex I refers to a sampling of low to medium frequencies and subindex h to a high 
frequencies one. Solid line would imply a = /3. See details in text. < 1. 

In order to avoid that these characteristics of the |' 0 / 3 | functions could make 
seem them belonging to a irregular family of functions ruled by a parameter, we 
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shall change the scale of the abscissae axis and we will draw our attention to the 
first and second derivatives of such functions in the new scale. The choice for the 
axis is the decimal logarithmic scale of variable a;/27r and the first derivative 
to look at is (logio |V’/ 3 |)'= Obviously (logip l^/’/sD" represents the 

second derivative 

(i(logio 

With this depiction it is clear that each (logj^Q \'ij^p\)' looks like a step whose 
height is greater as /3 approaches 1 from 0, evolving the figure from a slope 
to a sharp jump around the same interval of frequencies, G (10“^, 10^), and 
being almost constant outside it. On the other hand the second derivative 
(logj^o is ^ smooth peak, almost symmetrical, and of contracting half¬ 

width, which maximum moves from Vmax ~ 10“^ to i^max ~ 2 x 10“^ as f3 goes 
to value 1.00 from value 0.30, (see figure and b). 

That is, there is a plateau for each function \'4>0\, at low frequencies, and 
suddenly it bends in a neighborhood of j/ ~ 10“^, till it reaches a constant 
potential decaying tail, which is reflected in the graphic of the first derivative 
as a very quick approach to its horizontal asymptote of value /3. The greater 
the value of /3 the quicker the approximation to the horizontal, i.e. for small 
values of /3 a larger sampling of frequencies is necessary to fit the function |' 0 / 3 | 
with a potential decay and to get a ~ fi. To quantify how long this 

frequency interval should be, we have split the domain of modulus of tjjp in 
two regions, one corresponding to the plateau and fall, which we discarded, and 
the other corresponding to the tail, being the mark at « 20 and the end at 
v = 500.0005 -denoted with subindex Z -, or at z/ = 10^^ -denoted h -, and we 
proceeded to adjust the tail to the generic function with M and a the 

parameters of the fit. (See figure]^. It is shown that both series of exponents, 
ai and a^, follow closely the line corresponding to an ideal exponent (3. The 
first series, that considering frequencies till v = 500.0005, shows a slight bulge 
in the interval 0.10 < /3 < 0.30 where the absolute error is greater although 
the relative error (/3 — ai)/f3 is a monotonously decreasing result in the path 
from /3 = 0.02 to /3 = 1.00. As we pointed out, this is a consequence of the 
slower relaxation of the logarithmic slope, (logj^Q towards its asymptote 

for small values of beta and it doesn’t indicate at all a qualitative change in the 
behaviour of the modulus of fy/gj that could break the trend of tail decaying as 
. Nor it suggests the need of a functional description for the low frequency 
part of \tjjp\ for small betas different than for bigger ones. The second series 
corroborates this though as the absolute error |/3 — a/ij is barely perceptible for 
the values of /3 = 0.02, 0.05, 0.08 or O.xO with x S {1,2,...,9} and j3 = 1.00. 
Needless to say that the relative error is even smaller than the previous case for 
each P, (see figure [^). 

Following this a question raises when an adjustment to a probe function, (w. 
gr. a sum Havriliak-Negami functions), is tried. Is the test function ’flexible 
enough’ to describe the whole range of frequencies? Or is it better to describe 
the original data with two ’maps’, one for low frequencies and the other for high 
ones and matching them conveniently in the frontier? 

The question is not trivial, on one hand because the weight of the low fre- 
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Figure 2: As in figure a) Logarithmic derivative of logj^Q vs. logjQ(ij/27r). b) 

Logarithmic second derivative of logj^Q |i/)^(a;)| vs. logig(aj/27r). c) Exponent a of tail for 
fitting \ipp(u])\ ~ Subindexes I and h refer to two different sampling ranges of medium 

and high frequencies, respectively. Solid line, again, signalizes ideal asymptotic tail a = fi. 
See text for explanations. /3 > 1. 


quency part competes with the weight of the high frequency part and the former 
distorts the potential behavionr we have jnst described, moving away from —/? 
the value of the exponent of tails. And in the other hand nnmerical samplings 
and manipulations, as in FFT, broke the tail trend and ronnded it off affecting 
conseqnently the obtained values of parameters. These issues will be addressed 
in the following pages by means of snms of Havriliak-Negami fnnctions fitted to 
data calcnlated almost symbolically. 

However some distinctions shonld be pointed out with the case 1 < /3 < 2, 
before to continue. Although both are the only ones with a positive density 
fnnction of times of relaxation for all the stretched-squeezed exponentials ElE], 
exp(—f^), (j.e. with /3 g( 0, oo)), they behave differently in the derivative and 
both show an obvions parametric discontinnity in the t—space. While for /3 < 1 
the minus derivative, , shows a singularity at t = 0, the case /3 > 1 

presents a regular and continuons zero valne though it pays the price of losing 
monotone decreasing behaviour of the former. Now it displays an increasing 
monotone growth and subsequently a decreasing monotone tail which saves the 
character of the positivity of the density function mi. They conld be cases 
of a very different appearance or even to show a break in the continuity after 
doing the Fourier Transform for both but the real fact is that the parametric 
family of the first derivative transform, ipp(uj), does not present any oddity. 
Nonetheless the instance /? > 1 shows, trnth is, a qualitatively different tail- 
to-body (j.e. tail to plateau and fall) response. While for the case /3 < 1 the 
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tail have the look of a natural extension of the fall, in the case /3 > 1 there 
is an abrupt change of direction between the body fall at medium frequencies 
and the high frequency tail. We draw on to the representations of (logj^Q \4’p\)' 
and (logj^g |'0/3|)" to explain the changes. (See figure]^. The first logarithmic 
derivative it is no longer a simple step, after the drop a narrow well is formed, 
(with the bottom at z/ < 0.8), and following this the line comes up and stabilizes 
at the selected value of —/3, the same of the asymptote predicted by the series. 
What it is remarkable since the functional series from we obtained the expression 
for potential tail it is not converging in this case. (See figure]^). 

As a direct consequence of this new geometrical characteristic of the first 
derivative, the second logarithmic derivative presents two peaks, one negative 
and the other positive. The negative is a parametric continuation of the family 
shown for the case /3 < 1, and has as a new property certain degree of skewness 
due to the presence of the second peak. Again the position of the maximum in 
OX axis grows in frequency and approaches to Vrnax ^ 1- It never reaches that 
’high’ value for the interval 1 < /3 < 2, again because the second peak. The set 
of second peaks, looked at as a parametric family, accrues their maxima near to 
Vmax 1 as > 2. They are smaller than their corresponding partners. (See 
figure]^ there in the boxes of the first and second derivatives of log^g IV’/sIj for 
sake of clarity, are not depicted all the curves available but we have selected 
the last one of the previous group, /3 = 1, the curves corresponding to j3 = l.xO 
with X S{1, 2,..., 9} and /3 = 2.00). The first peak corresponds to the fall of the 
body of the function log^g \'>Pp\ and this one is even steeper than in the previous 
curves < 1 as the depth of the minimum grows in magnitude. The second 
peak, the positive, smaller in magnitude than the preceding negative, is a new 
phenomenology not seen in the previous case. It indicates a fold upwards in the 
curve logj^Q in a very narrow region of frequencies as the short half-width 
of the peaks is developed in an interval with upper limit not beyond ^ 2 as 
/3 —>■ 2. Of course for values of /3 < 1.40 the half-width is greater than this limit 
value but the height of the peak is so small compared to greater values of /3 that 
the sidestep of the curve logj^g in this region is barely perceptible. (See figure 
Beyond this range of frequencies {i.e. v > 10) the ’activity’ of the second 
derivative, (log^g |' 0 / 3 |)^^ is negligible and the function is mostly of a potential 
nature. In the figure [^, as in the former case a graph of the exponent of tail vs. 
P parameter is shown. We hold the same notation as before, ai for the result 
of the adjustment of Mv~°‘ to log^g in the interval v € (20.0,500.0005] 
and ah for the fit in the range of z/ S (20.0,10^]. For the naked eye is quite 
difficult to distinguish the trend of these two sets of parameters from the result 
a = P, (solid line in the figure), what it points to an early high frequency tail 
very different in behaviour from the main body of the spectrum. 



Figure 3: Left panel, a)-. Parameters a and 7 got from adjustment ^ A* HNa,'y{uj). 

The coarser step sampling, and A = 1, is for mesh ri. Finer steps are those of r 2 , {A = 1), 
and r 2 b? ^ 7 ^ !)• Right panel, b)\ Products a. • 7 for the previous parameters. Dashed lines 
would correspond to situations a. • 7 = (3 and = 


2 Common approaches 

2.1 The asymptotic trends of 1-HN approximation 

2.1.1 The stretched instance /3 < 1 

We have drawn two conclusions from the asymptotic behaviour of \il}p\. Firstly 
from a point of view of reconstruction in t—space only an small interval of low 
frequencies are needed to recover the main traits of the relaxation. Approxi¬ 
mately the most significant part, -plateau, bend and fall-, happens in the range 
of zy G [0,1] for most the betas and the tail almost stabilizes to its asymptotic 
functional form before v < 10. And second this tail performs quite well as a po¬ 
tential decay of exponent (3, i.e. \ipp \ « Mv~^. Consequently, keeping in mind 
this two important features, and with the intention to design the best possi¬ 
ble approximation, we shall test how the Havriliak-Negami function it sticks to 
the data, IV'/?!) when we optimize the error between them. In short we wonder 
how the approximation in area distorts the desirable common path of the tails 
and reciprocally how the high frequency concomitance of the tails spoils the 
very low frequency description of a function, which is steadier than its 

Havriliak-Negami alternate. 

To this end we made three related adjustments of the data in the interval 
V G [0, 500.0005] to one Havriliak-Negami function. In the first case the sampling 
step in frequency interval was 6v = 500/999.999, that is a purged sample of 
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the points used in the second and third cases, both with regular stride Sv = 
1/999.999. Thus the influence of the body is practically eliminated and the tail 
is the responsible of the values of the parameters obtained for each beta. In 
the second case a standard fit is done, and in the third one the initial condition 
HNa.-y{uj = 0) = 1 is relaxed to a free value, i.e. ~ HNa^~f{uj) with A 

to be determined in the optimization too. The idea behind this is to allow the 
Havriliak-Negami approximant compensates its lack of suitability in the very 
low frequencies, (there it doesn’t measure up, given its modulus is below I'/’/sDj 
and to soften the influence of the tail in relationship to the body part of data. 
The result, for /3 < 1, points to the sensibility of a and 7 parameters to small 
changes in the approximation of \il)p\ area by A* HNa^-y(ui), -even with so small 
variations of A as a maximum 6% around /3 = 0.32-, and the sturdiness of tail 
influence on the asymptotic trend of HNa^j. 

In figure to the left, we depict parameters a and 7 for the three referred 
cases, and in the right box the corresponding products a-j are graphed. In both 
boxes the subscripts ri are for the purged case, the r 2 ’s are for the common fit, 
ipp « HNa^ry, and the modified case il)p A* HNa,-y is denoted with r2b- 

We can see how ari{j3) and 7ri(/3) have a monotonous increasing behaviour 
and the curves only touch themselves at /? = 1. A changing picture when a finer 
mesh is used and the very low frequencies are added to the data, now the curves 
ctr 2 {P) and 7 r 2 (/ 3 )) still increasing ones that meet at /3 = 1, cross one to the other 
within the interval (approx, at 0.32 < /3 < 0.35). And the same happens to the 
extended approximation A*HNa^j, though the crossover of ar2biP) and jr 2 b{/ 3 ) 
it is shifted to the right of the former (0.42 < j5 < 0.45). The ar?{l3) curves don’t 
change qualitatively, as their all first derivatives keep monotonously decreasing, 
but the 7 r?(; 0 )’s do inasmuch as the first derivative suffers a transition from 
strictly decreasing function, (ri), to one with a minimum (i.e. it decreases, 
halts and then increases for cases r2 and r2b)- Besides for both families, a and 
7 , an apparent transition from ri to r 2 is observed in case r2b- It suggests 
an enhancement of the weight of tail in detriment of plateau and fall weights, 
what it sounds logic insofar as the extended condition frees the constraint of 
approximating the low frequency zone by the correspondent region of Havriliak- 
Negami function, and this allows a greater influence on the remaining part by 
the data of high frequency. That is to say we do not provide to HNa^-y{oj) with 
a new functional flexibility, this should be achieved for example with additional 
terms, by contrast now the fit is obtained by coating, instead of adapting, the 
plateau of data. All this sets free the function tail to adapt itself to respective 
data tail. Consequently the appearance of transition that shows r 2 b, stresses 
how carefully the design of the approximant, the method of optimization, the 
data absent of errors, and of course the frequency domain, should be done or 
chosen. 

This is also seen in the graph of the tree products a ■ 7r'?(/3), (figure right 
panel). The first one, ol ■ 7 ^ 1 , follows heavily a straight line y = (3, indeed a 
regression to y = 13°' gives a = 1.035 with a correlation of c.c. = 0.999988. The 
second and third cases (r 2 and r2b), give a = 1.143 with c.c. = 0.998633 and 
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Figure 4: Parameters, ri ^2 arid A in the panel a, and in the h one a. 1^2 and 71 , 2 ? for an 
approximant as in Eq. ^ obtained adjusting the function to, (ri), sampled 

data of in the range of medium frequencies. /3 < 1. Observe the similarities with a fit 
of same kind for high frequencies. Only A (medium) shows a behaviour more akin to low 
frequencies than to high ones. See figuresandDotted lines are just an aid for the eye. 


a = 1.113 with c.c. = 0.999004 respectively^ In all three cases the dependency 
deviates from proposed law y = (3 1-23 miiniiii], it is not much to say about ri 
since it holds data almost exclusively from tails, nevertheless r 2 and r 2 b main¬ 
tain, by sampling, the information needed to reconstruct the original function to 
large times (i.e. low frequencies), and is quite significant that the approaching 
to the mentioned law is different according to the value /3. When 0 < /3 < 0.6 
there is no closeness at all, and if 0.6 < /3 < 1.0 the agreement is only due to 
the similar functional dependence, /3“. The point here is, as /3 evolves from 0 
to 1 the (log^o \'4’^\y function takes form as an abrupt step which implies for 
1 ^, 31 a greater width for the plateau and a lesser magnitude and role for the tail, 
therefore while adjusting ipp to a Havriliak-Negami function, a progression from 
tail dominated, {y ^ m * f3, m < 1), to body dominated fit, {y ^ /3“, a > 1), 
will take place in product Logically this transition will be spoiled if the 

tails for small beta values does not develop properly and they are rounded off at 
the limits of an w—space finite interval, such thing happens for the FFT and is 
expected for an algorithm that erodes sharp characteristics whilst it determines 
the density of relaxation times |53L [SHI IHO] • In such circumstances each different 
procedure used will change the beta dependence of a • 7 in a particular and own 
way gaEsniniizi]. 

^As a matter of fact a better expression, 1 — (1 — /3)°', is possible though it would make 
more difficult the comparison to formulas like found in the scientific literature. 
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2.1.2 The squeezed instance /3 > 1 

The linear approach to a ■ 7(/3) nevertheless fails, for same kind of samplings 
Ti, r 2 and r 2 b, when 1 < /3 < 2. In this region of control parameter /3 a 
good empirical fitting will he y = /3+Ai(/3 — l)^-^+A2(/3 — +A'i{j3 — 1)^° 

with a correlation of 0.999999 for the three cases. Obviously the magnitude of 
parameters will vary being the look of ri different of r 2 or r 2 b- While at first 
sight ri could seem a deviation from linearity easily corrected with a simple 
power of (/3 — 1) -it is possible to correct j/ ~ /3 with a quadratic power to 
obtain a good agreement- the remaining two cases need more higher power 
terms to give the impression of a rapidly growing function. Besides in all cases 
Ai^ri <<Ai^r 2 b < with i = 1,2,3 and A 2 bigger than the other two 

A' s, allowing the interpretation of r 2 & as an intermediate case between a tail 
dominant adjustment and a body dominant one. Exactly as in interval 0 < 
/3 < 1, though with a great difference, the difficulty to obtain a behaviour that 
resembles the asymptotic trend of data y ^ (3. 

This hindrance should be borne in mind when making an approximation to 
tjj/}- As there is a sudden change of direction in the modulus from the dropping 
to the tail around « 1 and a very fast decreasing of magnitude, any function 
showing monotony in decreasing and lack of oscillations or ripples will have large 
difficulties to follow the data from low to high frequencies. Therefore our choice, 
a description by a simple Havriliak-Negami approximant, won’t keep on track 
with tail inasmuch as the low frequency part of data is heavier and prevailing 
than the high frequency part and the fitting function too mild for dodging the 
twist presented by the second logarithmic derivative. Consequently, even when 
a residual quantity of data of such frequencies are present in the sample, the 
optimization process will balance towards a selection of parameters describing 
mainly the plateau and rise height of step rather than the tail. This is what 
happens with case ri which deviates from linearity y = (3 and of course explains, 
when the whole sample of data is used in r 2 and r 2 b, the explosion of product 
a ■ 7(/3), for 1 < /I < 2. 

2.2 The asymptotic trends of 2-HN approximation 

While the approximation of tpis by means of one Havriliak-Negami function de¬ 
scribes properly tail decaying of data for /3 < 1, -although the asymptotic ten¬ 
dency of approximant suffers distortion because the matching to plateau zone-, 
it is clear from scientific literature that it is not enough to give an account of 
shape for the low frequencies uniMiiTniiHT]. Neither is the case for /3 > 1 as the 
advantage of picturing the high frequency asymptote is lost because tail weight 
in the optimization process can not offset plateau and drop predominance. How¬ 
ever functional proximity can be improved adding a new H-N term for adapting 
to EFT of finite size temporal data [371 Eg. Is our intention to show that the 
same is true for the numerically calculated complex numbers with an 

added feature, a right portrayal of the tail. Something that was not previously 
attained inasmuch as a finite window in the EFT distorted by convolution the 
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Figure 5: As in figureparameters {ai^2i 7l,2i'ri^2i A} are calculated for a double approx- 
imant of Havriliak-Negami type. See Eq. The frequency of data employed reaches up to 
V ~ 10^^ with /3 < 1. Lines are mathematical functions designed to describe these curves as 
dependent with variable /9. See tables and in appendix. 


high frequency part of spectrum. Our goal is then to show how this two terms 
approximant gives a good picture of both low and high frequencies. And how 
each term shares its contribution to the global approximation, a pertinent ques¬ 
tion because the roles of both terms will not be equal for cases /3 < 1 and /3 > 1. 


2.2.1 The stretched case /3 < 1 


We have used the series of data till v = 500.0005 in both versions r-i {i.e. 5vr2 = 
1/999.999), and ri {Siyri = 50Q*5i'r2) and the logarithmically homogeneous one 
fsi (Sui^si = 9 * 10“““^), until ly = 10^^, or v = 10^, for ,5 < 1 or ,5 > 1, 
respectively. The two terms approximant already described which will adjust 
data is 


2 

Ap2-ff^^a,7.r,A(w) = 

S = 1 


A. 

(1 -I- (iTsUj)^‘P‘’ 


( 1 ) 


with Ai = A and Ai -f A 2 = 1. 


Mesh sampling (ri) 

In figure are presented the fitting parameters for case ri which is dominated 
by tail values in medium frequencies as the low frequencies are mainly absent 
by decimation. The right panel describes the responsible ones for the potential 
behaviour of tails, ai ^2 and 71 , 2 , and the left panel describes the characteristic 
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Figure 6 : Parameters {oi, 2 i 7l,2i 'ri^ 2 i A} obtained after minimizing ||V'/ 3 | ~ 

|■4p2^^^^cI,'Y,T,A(<^)lI in the interval of frequencies v S [0,500.0005], with hv = 1/999.999. 
/3 < 1. Solid lines are analytical rough estimates to describe the parameters. See tables]^ and 
l^in appendix. 


times of both functions, ti. 2 , plus share coefficient A. In the left panel is conve¬ 
nient to take notice of a slowing varying lambda in interval 0.02 < /3 0.6 with 

a value near to 0.5 which suddenly drops to zero from /3 « 0.60 to 1.00, also 
is remarkable the difference in magnitude between ti and T 2 . While never 
exceeds the value of 1 . 1 , ti requires of a decimal logarithmic scale to be shown 
in the same graph. It gives the idea of a dominant function associated with 
T 2 only corrected by the one associated with in the necessary zone of low 
frequencies, an oversimplified picture because is more accurate with values of /3 
near to 1.00 than with values near to 0.02 as both functions decay very slowly 
in this range of parameter /3, due to small values of a\ and a^. However it is a 
good frame of reference to talk about the contribution of each i 7 /Vs=i, 2 (w), so 
we’ll refer recurrently in these terms to them. 

In the right panel two features should be noticed, first and 72 do not cross 
one each other in line with what we saw while using only one Havriliak-Negami 
function to describe pruned data. And second in consonance with the fact that 
A diminishes in value approximating to zero and the numerical optimization is 
locked by only one HN{uj) function, ai and 71 both present an apparent discon¬ 
tinuity at ,5 = 1. Although not all of these two functions is wrong as they cross 
themselves a property of a and 7 already shown in the case of approximating 
data, (low frequencies included), with only one Havriliak-Negami function. This 
fact enforces the idea of considering HNi{uj) a backup for HN 2 {uj) while fitting 
functions to data in low frequency zone. (See figure . 
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Mesh sampling (rgi) 

With the aim of testing this hypothesis, namely the residuary influence at low 
frequency of ancillary Havriliak-Negami function {HNi{uj)), a further adjust¬ 
ment to data with larger interval of frequencies is done. Now the sampling called 
logarithmically homogeneous, rsi, is used as previous data have a very fine step 
Sv, which makes computationally hard any optimization in a sensible amount 
of time. The interval of frequencies is now v G [1.00,10^^] for /3 < 1 with the 
aforementioned step Svi^si which changes in each decade of frequency. After cal¬ 
culations a new layout springs out as result for the parameters of approximant, 
this is shown in figure We observe important differences, comparing to latter 
result, for every parameter. For share coefficient lambda the quick decreasing 
starts sooner, nearly /3 ~ 0.4 and has opposite curvature to the prior case, be¬ 
sides the values attained in the interval 0 < f3 < 0.4 are comprised between 0.5 
and 0.7 and do not seem steady as before. The curve for T 2 also suffers mod¬ 
ifications, now is smoother and starts near zero for /3 ~ 0 and reaches one for 
j3 = 1 . This is in consonance with the errors and the final optimization for r 2 
sampling case and suggests that the former curve for ri suffers from an strong 
source of error, possibly an interference between HNi{uj) and HN 2 {uj) in the 
initial interval 0.02 < f3 < 0.18. A hint to explain such situation is to realize 
the opposite results of ti and ai. While Ti would imply a very quick decay of 
HNiiuj), the small value of ai would slow down such decay, and would yield a 
non negligible contribution, something not to be expected from greater values 
of ai- 

This balance allows a slight rivalry with the values of HN 2 {uj) also dominated 
by a small a 2 - As a consequence of this equilibrium large fluctuations of T 2 in 
the neighborhood of ~ 0 make sense. Finally we should add that interferences 
in the vicinity of /3 ~ 1 between HNi{uj) and HN 2 {uj) also occur affecting ri, 
because for this curve a discontinuity happens at /3 = 1 since ri(/3 = 1) = 1 and 
ri(/3 <!)<!. And ti is not the only parameter with such discontinuity, ai and 
7 i also take values quite far from their ideal amount of 1 near /3 ^ 1 . (See figure 
[^. This pathology, partially consequence of a negligible A, points to a minor 
and complementary role of HNi{uj) describing ^pp<Cl{uJ) tail. A trend which 
is the opposite when the beta values are small, as the change of almost every 
parameter curve in the neighborhood of /3 ~ 0.4 seems to suggest. It is worthy 
to note the coincidence of this turnaround with the fact the tails are not fully 
developed at very high frequencies for lesser values of beta than 0.4, at least 
when a comparison with greater values, (v. gr. /3 > 0.6), is done. See figure 
graph a. To finish we should like to remark the quasi constant behaviour of 
71^2 in the interval 0.40 ^ /3 < 1 and the quasilinear one of ai _2 in the whole 
interval of beta, making the relaxation of tails for HN 2 {uj) nearer to a Cole-Cole 
type, (7 = 1), than to a Cole-Davidson one, (a = 1). By contrast relaxation for 
HNi{uj) still remains Havriliak-Negami in both cases ri and Vgi though they are 
very different as a family of curves, (compare 71 in figures and [^ , inasmuch 
as 7 i in ri is an increasing function for all the beta domain and 71 in rsi holds 
itself as quasi constant in the mentioned subinterval /3 G (0.4,1.0]. Surely the 
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lesser weight of medium and low frequencies in r^i with respect to ri takes its 
toll. 

Mesh sampling (r 2 ) 

Finally an adjustment to an r 2 sampling which includes many points corre¬ 
sponding to low frequencies is done, the result is that the roles of character¬ 
istics times Ti and T 2 are interchanged. Now ti scale does not surpass the 
value of one. For betas near to one ti is around 0.7 and at /3 = 1 is dis¬ 
continuous, besides after a maximum located at /3 « 0.62 it decreases to 0.1 
as /3 —)■ 0.02. Meanwhile T 2 is a monotonous decreasing curve from /3 ^ 0 
to (3 = 1.00 which can not be represented in the same scale as before unless 
decimal logarithms are taken and rescaling is done. (See left panel in figure 
. This suggests for the whole range of beta two scenarios, one in which the 
first term of the approximant, HNi{uj), losses progressively importance as beta 
grows. It starts at ,5 « 0.62 and ends when the approximant becomes degen¬ 
erate Ap 2 HNa^^^r,\^o{^) « -f^^{a 2 . 72 ,T 2 }-s-{i,i.i}(‘*^). at /3 = 1.00, (j.e. tends 
to a Debye relaxation as A —?► 0). This feature avoids the proper numerical 
determination of all parameters of HNi{uj) and explain why ti is discontin¬ 
uous, the optimization is locked by the suitability of HN 2 {oj) and no further 
contribution of HNi{uj) is possible. On the contrary, in the second scenario, 
(0.02 < f3 < 0.62), lambda, which takes values near to 0.4, shares out equal 
amount of contribution to each Havriliak-Negami term. Now the difference in 
contribution between both terms is due to the very distinct scale of time for 
Ti and T 2 , making one to be the complement of the other in this fit, while at 
the same time the tail asymptotic behaviour is mostly preserved. This situation 
allows to the approximant Ap 2 H^ description for a wide range of 
frequencies without excessive gap with data, (See figure]^. Pairs ai, 7 i 

and 02,72 show a similar behaviour, -for each pair, the curves cross themselves 
and have similar shapes to the case of one Havriliak-Negami approximant (see 
figure [^-, what is also consistent with a detailed description of plateau and 
drop, {i.e. the low frequencies), by jointly both Havriliak-Negami functions. 
(See right panel in figure . 

The empirical functions for the adjusting parameters ai, 2 i 7 i, 2 i ti ,2 and A 
are given in tablesandfor long frequency tails (vsi), and in tables[^and|^for 
low to medium frequencies (r 2 ), all of them are written for interval 0 < /? < 1 
of control parameter. 

2.2.2 The squeezed case (3 > 1 

However a different kind of results are obtained for the interval 1 < j3 < 2 . We 
already pointed out how a sudden change in direction makes bend upwards the 
curve |'!/'^(a;)|, which stabilizes in a potential tail at medium frequencies, after 
a very sharp dropping in the low to medium range of them. As a result the 
transition from plateau and rise height to tail is marked clearly by the decimal 
logarithmic curvature -which presents a chasm and a peak-, and consequently 


16 




Figure 7: Parametric curves { 01 , 2 , 7 i, 2 ,'ri^ 2 , A} of a two-terms approximant, (see Eq. [^, 
for data, Iti’/sl, picked from interval v £ [0,500.0005] with 5v = 1/999.999. /3 > 1. Solid lines 
are mathematical expressions tailored to adjust the curves. See tables and in appendix. 


the adjustment to data by a Havriliak-Negami approximant is made more dif¬ 
ficult, as the latter does not show such an extreme changes in curvature. The 
weight of low frequencies is bigger than the data corresponding to tail because 
their respective magnitude and total number of points, (for sampling r2), thus 
the parameters obtained after optimization are those of a quick decaying func¬ 
tion not following the natural asymptote]^ (See figure]^. Accordingly a global 
description of data would require necessarily a different two terms approximant 
of Havriliak-Negami type for the tail and every /3. To obtain this we would 
proceed as before with an additional fit to a sampling like ri or r^i as both 
are pruned of low frequencies. In these cases they are sets of data weighting 
mostly the medium to high frequencies and the corresponding result will de¬ 
scribe the tail asymptotic behaviour properly at the expense of not describing 
at all the plateau and its environs. Nevertheless the analysis 01,2, 71,2, '^1^2 and 
A as empirical functions of /3 deals with remnant information of this neglected 
important part of the spectrum. The triplet {02, 721 T2} is in both samplings, 
ri and describable by means of smooth functions {v.gr. combination of 
exponentials and polynomials) and also is the case, except for some occasional 
dimple in for the share coefficient A. (See figure]^. Nonetheless it is not 
such a way for the triplet {ai, 71, ri}. Though it could be possible to adjust ti 
with just one mathematical expression and the errors with data points wouldn’t 

®Due to the commented features of the function \ipp\ near i' ~ 0.9 the influence of clean 
data in the adjustment to a Havriliak-Negami approximant it is not very different of same fit 
to spoiled, (by EFT), data. In consequence the shape of empirical parameter curves do not 
change significantly in either, theoretical or EFT, case. See Figure [t] 
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Figure 8: Parametric curves { 0 : 1 , 2 ? 7 i, 2 » ■ 7 "i, 2 ? of ^p 2 -^^^a, 7 ,T,A function approximating 
data up to ~ 10^, in a homogeneous-like logarithmic pace. (See Eq. ^and text). Notice 
the big dispersion of experimental points 01 and 71 , and the correlated corresponding error 
of 02 and 72 . It is remarkable that, even with this scattering of the parameters, the products 
Ql • 71 and 02 ■ 72 hold themselves extremely attached to a linear trend. (See figure]^ and 
d). Solid lines are rough proposals based on constraints given by same products. See tables]^ 
and in appendix. 
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be large, it won’t be the case for ai and 71. (See right panel in figure]^. These 
last two curves are crooked, or even broken, in many points and do not show a 
defined trend, besides this happens for both cases ri and but each in their 
own ways. Thus only a piecewise description is admissible for them. However 
this is not the treatment we are looking for the parameters as it would require 
an interpolating table, and certainly many errors, instead a simple formula. But 
not all is lost, paradoxically the product oi • 71 is well behaved and can be fairly 
approximated to oi • 71 ~ 2/3 in case and to oi • 71 ~ — 0.22 + 2 . 06/3 in case 
ri. And meanwhile the product o.’i ■ 72 is even closer and more correlated to a 
straight line, (he. 02 • 72 — /3 with very small errors in both cases), than the 
previous product. The evident conclusion is that in spite of the bad conditioning 
and many errors finding (X\ and 71 out the asymptotic tendency of tails is well 
described with a two-term Havriliak-Negami sum. On one side the second term, 
HN2{uj), dictates the tail behaviour which at large ly coincides with the trend of 
real data, and on the other hand the first term, HNi{uj), -quickly fading-, gives 
a proper placement of the curve near such data. And for this last constituent, 
notwithstanding the poorly defined trends of a— and 7—type parameters, the 
stability of constraint oi • 71 ~ 2/3 allows us to put forward an ’ansatz’ for their 
dependence on / 3 , which finally enables a finer classification of high frequency 
spectral data. This will namely: 01,2 ~ 0 {/ 3 ), 71 ^ 0 ( 2 ) and 72 ~ 0 [X). 

In Figure]^ all a -7 products already commented for 0 < /3 < 2 are depicted. 
Subscripts numbers 1 and 2 stand for the first and second HN-f^uj) term in 
each approximant, and letters h and t describe the adjustment to r2 and r^i 
samplings, respectivelyH 

Also in tables EE and mathematical adjustments to parameters oi 2? 
7 i,2i ti 2 and A are recorded for 1 < /3 < 2 . In table and table the formulae 
represent a quite close match to the empirical points which define the approx¬ 
imant for low frequencies, the counterpart for high frequencies is in tables 
and[^ However in these last two tables the intended formulas for point curves 
ai, 7i and ti are just rough estimates as such curves show breaks and outliers. 
The remaining point graphs, (he. a 2 , 72) t' 2 and A), exhibit nevertheless an 
appropriate fit and as a consequence the mathematical expressions which sub¬ 
stitute them give a good description of the high frequency behaviour (he. the 
asymptotes) although possibly with a misplacement due to errors describing the 
set of parameters {oi, 71, ti}. 

3 Discussion 

By comparing the parametric curves {ai,2) 7i,2) 'ri,2) A}, (see figure]^, with those 
of a previous approximation to data '0,9 (w) which undergo convolution with 
a sine function, one can see a strong distortion in curves ai^2 nnd 71 2 for 
low values of /3 < 1 . (See figures 3 and 4 in reference [ 53 ]). In that case 
they didn’t tend to zero as it happens in our present situation for 017, nor 

^Denoting h the head, i.e. plateau and drop, and t the high frequency tail, i.e. potential 
asymptotic trend. 
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Figure 9: Products 01-71 and 02 • 72 of the asymptotic parameters, gotten by optimizing 
Ap 2 HNa ry ,,- \ ~ -i/)^ in samplings r 2 , —subindex h: head—, and r^i, -subindex t: tail—, along 
the whole interval of variable 0 < /5 < 2. The solid line inside boxes a, b and d represents 
a ■ 'y = /3, while it is for o ■ 7 = 2/3 in box c. In boxes a and b is also depicted the tendency 
Q ■ 7 = ^t.23 ^iiii a, dashed line. 


retained a value below 1 for 71 , whithin an interval of small values of beta 
{P < 0.25). And of course neither of them allowed us to consider constraints 
of type (a • 7 ) 1,2 ~ C>(/3) valid ones. However a successful reconstruction of 
the Kohlrausch function was still possible by means of inverse FFT for small 
(/3 = 0.15) and large (/? = 0.85) values, despite these facts. 

Some circumstances played a role to compensate the until then unknown dis¬ 
tortion of beta-parametrized curves {a, 7 }i, 2 - Among them: the serious prob¬ 
lems in the medium-to-high frequencies interval, v S [10,500], were overcome 
easily with a proper damping of phase in the approximation, while the prob¬ 
lems in low frequency range -consequence of discrepancy between data and the 
approximant pair- were shadowed by the finite size of time window through the 
smearing of peaks. Also the tails of both functions differed so slightly in a finite 
interval that made imperceptible the disagreements in reconstruction at short 
time scale. 

Nevertheless what it was really good for the purpose of ignoring the distance 
in function space between the Fast Fourier Transform and the double Havriliak- 
Negami approximant, -i.e. its adaptability thanks to the many parameters-, 
played against a proper description of the theoretical Fourier Transform. The 
main three barriers cited which oppose to the perfect fit of HN approximant and 
FT data, will suffer dissimilar transformation through the process of convolution 
of data with a finite window, by the simple reason of being differently developed 
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head, drop and tail along the same frequency interval. So, the change of shape 
induced by convolution, not only in plateaux but also in the tails where the 
potential behavior uj~^ was rounded due to finite size of spectrum, will be 
different with the value of /3. 

In consequence the functional dependence of {(a, 7, r)i^2j ^}(/3) for the HN 
functions has to result distorted with respect to the desired result of same family 
while follows theoretical FT data. Such deviation will depend on j3 inasmuch 
as the distortion of the FFT is not equally distributed with this parametric 
variable. Obviously a lesser drop tail-to-head in the functions with small /3 
values wore a greater adaptation of the approximants to the convolved data, 
loosing some of the original characteristics of theoretical data. 

4 Conclusions 

It is common practice while studying complex systems to analyze their relax¬ 
ations in time as well as in frequency. Unfortunately there are not often at hand 
short and compact expressions corresponding simultaneously to the mathemat¬ 
ical formulation of a same phenomenon in both spaces. This article is focused 
towards the approximation of Fourier Transform of certain Weibull distributions 
by Havriliak-Negami functions, studying their fitting instabilities and suitabil¬ 
ity. As the suggested models possess many parameters this first part delves 
thoroughly into the functional dependence of them with shape parameter f3. 
Also we discuss the influence of data variation and extension on the resulting /? 
dependency, i. e. the sensibility of fit with data errors or distortions. Finally 
we will end this part establishing additional restrictions on parameter behaviour 
which are deduced from asymptotic trends of the theoretical functions. 

When the values of '0/3 (‘j-') ar® sampled properly as to balance the weight of 
low frequency data with the more numerous and tiny values of high frequency 
tails during a fit or optimization, and besides such data are not distorted or 
biased by finite time windows in the Fourier transform of the original function 
— ^ or other numerical artifacts as consequence of switching from time 

to frequency space |58l[70llMllM| . it will be true that a unique Havriliak-Negami 
function shall not follow other asymptotic condition a-y = /3^ different of d = 1. 

Also whether two Havriliak-Negami functions are used as approximant the 
asymptotic constraints will be (a ■ 7)1,2 = / 3 , for /3 < 1. While (a ■ 7)1 « 2/3 
and (a ■ 7)2 = /3, for 2 > /3 > 1. Nevertheless whenever the heads, (ly < 1), are 
involved and the tails, (1/ > 1), are not sufficiently important to talk about a 
fully developed asymptotic trend, (see figures and [^ , the constraints for a ■ 7 
won’t be fulfilled. In such cases (a-7)i,2 ~ 0 {l 3 ) for /3 < I, and (a-7)i,2 ~ 
with d > 2 for /3 > 1 . 

These constraints should clear up our understanding of the fitting imped¬ 
iments found at low frequency while the model is being built. They suggest 
for slight changes in the uniform model of an approximant with two terms to 
give account of the whole range of frequencies. Unfortunately the accuracy of 
such approximations is restricted to the description of modulus of the objec- 
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&A exp (X;; 


s{l - /5)=)exp(-M/3<i) 


^exp as(l-/3)^) exp(-M/3fi) 


Parameters 

Constants 

ai 

0-2 

7 i 

72 

A 

1.68125 

= 1.000 

1.6738 

= 1.000 

M 

6.00482 

2.14961 

4.95993 

2.49211 

d 

0.712224 

0.580585 

1.17881 

0.651276 

ai 

-30.8645 

0.56705 

-5.87091 

0.585986 

Cl2 

75.8761 

-2.83736 

13.4251 

-2.92872 

as 

-44.5841 

2.97408 

-6.94891 

3.0715 

Corr. 

0.999858 

0.999997 

0.997978 

0.999256 


Table 1: Approximant of tails, case /3 < 1: Formulas used to adjust {cn, 2 i 7 i, 2 } and their 
fitting constants 1881 . Corr. = correlation coefficient of nonlinear fit. For oi and 71 , /3=1 is 
an outlier. See right panel at figure]^ 


tive function as far as our primary purpose is to understand the modifications 
required over HN functions to fit data in a better way. An additional work to 
fully optimize our models in the complex domain should be done by means of 
the information obtained from the new family of relaxation functions we have 
introduced here. Work is in progress towards such direction. 


Appendix 

Here are given the tables with formulae and parameter values corresponding to 
mathematical adjustments of point data in figures to 

Case 0 < /3 < 1 

Tablesanddescribe a high frequency model and tables]^ anddo the same 
for a low frequency one. 

Case 2 > (3 > 1 

Now the first pair of tables, an di are those describing a two-function approx¬ 
imant in the domain of low frequencies, and tables and the corresponding 
to the high frequency model. 
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-0.65664 + . , . 



+ Es = l |^a(l - f>) + Bs{l - /3)^| exp(-Ma(30-2) 

/3^ + EI 1 iis(l - 13}“ 

Bxp(-M^‘i) Eb=i cs(1 - I3)“ 



Parameter 

Oonstants 


logic 


Al 

-22.4446 

Bi 

-6.64276 

Ml 

-0.233028 

^2 

8231.54 

B2 

-7799.77 

M2 

5.59666 

Corr. 

0.999984 


Parameter 

Constants 

r2 

d 

2.96596 

bi 

4.54749 

b2 

-32.4329 

b3 

552.577 

bi 

-3159.47 

h 

9513.52 

be 

-13492.1 

br 

7139.92 

Corr. 

0.998613 


Parameter 

Constants 

A 

M 

4.14933 

d 

3.1579 

Cl 

3.43714 

C2 

-12.6851 

C3 

35.8224 

C4 

-46.5211 

C5 

20.5172 


Corr. 0.998905 


Table 2: Approximant of tails, case /S < 1: Formulas used to adjust {ti, 2 }A} and their 
fitting constants. Corr. = correlation coefficient of nonlinear fit. See left panel at figure]^ 


/3 exp |^(Es^l “s(l - exp{M/3'^)J 


Formulae 









Parameters 

71 

72 

Parameters 

OL\ 

CX.2 

Constants 

Constants 

M 

4.54457 

3.4752 

M 

3.6903 

2.42224 

d 

1.62482 

1.47877 

d 

1.4618 

1.03791 

A 

0.999202 

1.00287 

ai 

0.0190435 

0.0962135 

ai 

0.0227205 

0.0401981 

0,2 

0.0794805 

0.193794 

-0.055399 

0.054449 

02 

oz 

0.572997 

0.761107 

0.79403 

1.29633 

0.3 

Corr. 

0.999950 

0.999942 



Corr. 

0.998473 

0.999654 




Table 3: Approximant of heads, case /3 < 1: Formulas employed to adjust {o:i, 2 » 7 i, 2 } and 
their fitting constants. Corr. = correlation coefficient. See figure]^ 



- 0)^ 

- 0) + 

- /3)2| exp(-Ms/30-2) 

exp[-M(l - ^)] E® = 1 cs(l - 6)“ 

Parameter 


Parameter 

logio 3-2 

Parameter 


Constants 


Constants 

Constants 


d 

0.785731 

Al 

42.4745 

M 

4.60773 

bo 

0.726519 

Bi 

-49.3268 

Cl 

4.80377 

bi 

0.160075 

Ml 

3.65441 

C2 

-12.6617 

b2 

7.16583 

^2 

66182.6 

C3 

78.3249 

bs 

-14.5815 

B 2 

-61430.8 

C4 

-120.494 

bi 

8.47342 

M 2 

12.7917 

C5 

89.8475 

Corr. 

0.997479 

Corr. 

0.999987 

Corr. 

0.996847 


Table 4: Approximant of heads, case /3 < 1: Formulas employed to adjust {ti, 2 , A} and their 
fitting constants. Corr. = correlation coefficient. See figure 
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Formulae 

l + = - 1)», 

e = 10-180. 

A + exp(MO - 1)1 

)E^=i - i)“ 




Parameters 

71 

72 

Parameters 

01 

02 

Constants 

Constants 

A 

1.90381 

= 1.000 

M 

0.116811 

0.165105 

M 

2.25005 

0.239754 

ao 

0.136577 

-0.0645705 

bi 

-4.09964 

0.948757 

ai 

-0.0269458 

0.0697225 

22.844 

0.634155 

&2 

0-2 

-0.00692586 

-0.0485374 

-29.9149 

0.692054 

63 

0-3 

-0.00572689 

0.0171122 

11.5881 

-0.484575 

64 

Corr. 

0.999249 

0.999926 



Corr. 

0.999956 

0.999993 




Table 5: Approximant for heads, case /3 > 1: Formulas used to adjust { 0 : 1 , 2 ? 71 , 2 } their 
fitting parameters. Corr. = correlation coefficient. The point at /3 = 1 is an outlier for a fit 
to 7 i and is removed. See figure]^ 


coV(/3 - 1) + = l - 1)^ 


exp[-M(/3 - l)]Es = i - 1)® 


Parameters 

Constants 

CO 

Cl 

C2 

C3 

C4 

Corr. 


- loglO !■! 

- loglO T2 

0.339472 

0.27042 

0.27034 

-0.1091 

0.0194173 

0.977286 

0.26079 

-0.884491 

-0.198826 

0.303991 

0.999853 

0.999961 


Parameter 


Constants 


M 

2.9177 

di 

2.80567 

d2 

-3.18276 

ds 

4.6885 

Corr. 

0.998497 


Table 6: Approximant for heads, case /3 > 1: Formulas used to adjust (ti, 2 ,A} and their 
fitting parameters. Corr. = correlation coefficient. See figure 


Formulae 

13(^1 + 

ai (/3 - l)j 



A 


(^l + ai{/ 3 -l)j 

Parameters 

Constants 

ai 

ai 

-0.127772 

02 

-0.00819886 

Parameters 

Constants 

A 

ai 

71 

=2 

-0.110735 

72 

=1 

-0.0102285 

Corr. 

0.981962 

0.994925 

Corr. 

0.833974 

0.096606! 


Table 7: Approximant for tails, case /3 > 1: Very rough trends decide formulas used to 
adjust { 01 , 2 , 71 , 2 } since noise is greater than signal at many points. The points at /3 = 1.50, 
1.52, 1.55, 1.80, 1.82 and 1.85 largely break the trend for a fit to oi and 71 consequently they 
are removed. See figure]^ Notice the strong linear feature of products (o • 7 ) 1,2 despite the 
noise (see figure]^. 
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, 




B + 

- 1 + - 1)" 

exp[-ArV^ - 1] - 1)" 

Parameters 

Oonstants 

A 

- logio ri 
-0.278737 

- loglO "^2 
=0 

Parameter 

Oonstants 

A 

bo 

1.69465 

0.165075 

M 

-0.504694 

bi 

-4.114 

0.219429 

Cl 

1.41581 

&2 

9.33908 

1.33322 

C2 

-1.49885 

&3 

-10.9942 

-2.39982 

C3 

0.572827 

&4 

4.63626 

1.20034 

Corr. 

0.998585 

Corr. 

0.996880 

0.999016 




Table 8: Approximant for tails, case /3 > 1: Formulas used to adjust {ti, 2 } A}, and their 
fitting parameters. The points at /3 = 1.12, 1.15, 1.78, 1.80, 1.82 and 1.85 are not taking into 
account while adjusting. See figure]^ 
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